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Abstract 

We derive semiclassical neutrino-electron transport equations in the collisionlcss (Vlasov) limit from 
the coupled Dirac equations, incorporating the charged and neutral weak current-current as well 
as electromagnetic interactions. A corresponding linear response theory is derived. In particular, 
we calculate the response functions for a variety of beam-plasma geometries, which are of interest 
in a supernova scenario. We apply this to the study of plasmons and to a new class of collective 
pharon resonance modes, which are characterized by u) < q. We find that the growth rates of the 
unstable modes correspond to a strongly temperature (oc T^T^) and linearly momentum dependent 
e- folding length of about 10 10 km under typical conditions for Type II supernovae. This appears to 
rule out such long-wavelength collective modes as an efficient means of depositing neutrino energy 
into the plasma sphere. 

PACS numbers: 05.60.-k, 13.15.+g, 14.60.Lm, 97.60.Bw 



1 Introduction 

Neutrino transport processes are known to play a major role in the energy-momentum flow powering 
the dynamics of Type II supernovae 0, B|. Generally, it has been that the collision-dominated 
aspects have been studied in detail, leading to substantial progress in the understanding of these 
stellar explosions, while, however, still leaving open some problems in the quantitative description of 
their spatio-temporal (hydrodynamic) evolution. 

An earlier work by Bethe ||], which introduced the idea of a modified in-medium neutrino dispersion 
relation and of a corresponding effective Hamiltonian, was used in the series of papers [|| to describe the 
collective interaction of an intense neutrino flux (from the supernova core) with an electron-positron 
plasma (the supernova atmosphere) of comparatively low temperature. Their tentative conclusion 
was that a particular induced plasma instability may be much more efficient than traditional collision 
dominated mechanisms, i.e., faster by many orders of magnitude, in depositing the neutrino energy 
into the plasma sphere || . 

However, the approach in these papers is subject to criticism, since there is no physical or formal 
justification (such as a hypothetical condensate) for the scalar 'bosonic' collective neutrino wave func- 
tion used. In particular, the implied quantum phase coherence of the neutrinos appears hard to justify. 
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Considering their incoherent thermal production and the effective duration or length of the 'beam', 
no bunching effects are to be expected. Moreover, it is somewhat hidden in their phaenomenological 
approach how the precise V-A tensor structure of the electroweak current-current interactions can be 
taken into account. This has also been pointed out in Ref. Q recently. Employing finite temperature 
field theory, Bento studies the excitation/damping of longitudinal electromagnetic plasmons (Lang- 
muir modes) in the electron plasma under the influence of a neutrino flux. His results indicate that 
this type of collective mode instabilities "... do not seem to be a viable mechanism of substantial 
energy transfer from neutrinos to a supernova plasma" |J. We confirm this. 

It is the purpose of our present work to systematically derive the transport equations for the 
neutrino-electron system from first principles (Section 2), as well as the relevant dispersion relations 
(Section 3). We introduce appropriate spinor Wigner functions, while deriving the detailed chiral 
structure of the neutrino Wigner function in the Appendix. Previously, only the phenomenological 
approach []|] or the perturbative finite-temperature field theory |], f?j were applied. Our general 
derivations may also prove useful for other astrophysical applications, such as those involving strong 
magnetic fields or, generally, neutrino transport under mean field conditions. 

In the collisionless regime, the results of Sections 2 and 3 allow us to investigate, in detail, the 
collective modes in the highly anisotropic neutrino 'beam' plus electron-positron plasma system. We 
find longitudinal and transverse plasmons, which are only perturbatively modified by the neutrino 
flux. 

Furthermore, we also find a new class of growing, as well as decaying collective oscillations, nonex- 
istent in isotropic equilibrium plasmas, which we name 'p/iarons'Q. They are caused by a resonance 
effect, generally at a frequency to less than the momentum q, due to the unbalanced neutrino momen- 
tum distribution, which is characterized by a finite opening angle with respect to the beam axis. We 
study such modes with the wave vector parallel to the beam direction (Type I pharons), as well as 
with the wave vector orthogonal to the beam direction (Type II pharons). 

Geometrically, the Type II situation corresponds most closely to the one where the two-stream 
instabilities would be expected to occur in ordinary plasmas. We investigate whether such instabilities 
are induced by the weak current-current interactions and, depending on the growth rates, may provide 
an essential contribution to the still partly elusive energy transfer mechanism in Type II supernovae 
(Section 3). 

The collective two-stream filamentation instability is well known to occur in ordinary plasmas due 
to the electromagnetic Lorentz force [§, ||. More recently, it has also been studied in the context of 
strong (color-electromagnetic) interactions, where two interpenetrating parton beams describe high- 
energy nuclear reactions, see [|l0| and further references therein. 

More generally, one may expect such 'hydrodynamic' instabilities in interacting many-body sys- 
tems, in particular, in plasmas with interactions mediated by the Standard Model gauge fields, when- 
ever the system consists of two or more components with considerably different momentum space 
distributions & |], [l(J. In these cases, perturbations which, loosely speaking, are transverse to a 
predominant collective flow, tend to be amplified by the collective feedback effect of the effective 
long-range forces of mean field type. 

Consequently, we are motivated in this study by the supernova geometry, where the radially 
outward streaming neutrinos interact with the electron-positron plasma, which may produce a variety 
of collective instabilities. 



After the island of Pharos, where the famous lighthouse of ancient Alexandria was constructed under the order of 
Ptolemeus II. 
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2 The Coupled Neutrino-Electron Transport Equations 

Our derivation of mean field transport equations will follow the successful strategy developed earlier 
for QED QCD [12], and hadronic matter [JD!]. The basic idea can be easily summarized as follows: 
Starting from the underlying field operator equations of the model under consideration, one converts 
these into corresponding Wigner operator equations, i.e. for the density operator in the Wigner 
representation. In the appropriate mean field approximation the latter can be converted into a closed 
set of Wigner function equations (cf. Section 2.1); furthermore, performing a consistent ^-expansion, 
the most relevant semiclassical (Vlasov type) transport equations for the coupled relativistic phase 
space distributions are obtained (cf. Section 2.2). - Presently, our notation and conventions follow 



those of Ref. JTl ] 



2.1 Prom Dirac's to Mean Field Quantum Transport Equations 

We add the effective local coupling terms, 

G F 



Ant - ^ 



^ H 7m(1-75)^ H ^ ( V(<v-CA-ys)^ (e) , (i 



to the free Lagrangian densities of the electrons and electron neutrinos (including their antiparticles) . 
They represent the weak charged and neutral current- current interactions of the Standard Model in 
the appropriate low-energy limit, with cy = \ + 2sin 2 #y[/ and ca = \, using standard notation ||. 
For [i- and r-neutrinos, with only the neutral current interaction contributing, cy,A cv,A — 1- Here 
the neutrinos are described by four-component spinors as well and we allow for a small but finite 
neutrino mass, taking the growing evidence into account [jlj]. Eventually, however, we will pass to the 
massless limit, since for our applications the masses of order eV are definitely negligible compared to 
the relevant energy scales of order MeV. 

The electromagnetic interaction of the electrons must be added to the interaction of Eq. (|l]). For 
the derivation of the transport equations, however, this does not introduce a new element. The 
corresponding modifications will be added at the end of this section, making use of the earlier QED 



results pi [. 



The resulting Dirac (Heisenberg operator) equations for the electrons and neutrinos, incorporating 
the interaction (|l]) in the mean field (Hartree) approximation, can be written in the form: 

| i7 . dx _ m (Q _ j(0 . 7 ( c « _ C W 75 )} ^(0 = o , (2) 

where I = e, v denotes the electron and neutrino case, respectively, and Cy A = cy,A, Cy\ = 1; the 
neutrino current has to be inserted into the electron equation and vice versa, as indicated by J^'- 1 
here. The V-A four-currents are defined by: 

jP^^(:^ l h^-c^ 5 )^:) , (3) 

where the expectation value of the normal-ordered product refers to the ensemble characterizing the 
state of the system, which will be specified in more detail later. 

Introducing the Wigner functions, i.e. (4x4)-matrices with respect to the spinor indices which 
depend on space-time and four-momentum coordinates [JO]: 

W$(x,p) = J ^^ e -^(: $f( x + y/2)^(x - y/2) :) , (4) 
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and with the ^-dependence made explicit, the currents of Eq. (||) can be expressed as: 



j(i) {x) = 9Z tv J d 4 p 7M (cg ) -c2 ) 7s)W«(x,p) , 



(5) 



with the trace refering to the spinor indices. 

Multiplying the Dirac equations (||) with the respective adjoint spinor and making use of Eqs. (||)- 
(|5|), they can be converted to the Wigner representation. This yields the coupled electron and neutrino 
quantum transport equations: 

(7 • K - m«) W®(x,p) = exp(-y^ • d p )J«'\x) ■ 7 (cj? - c { ^ 5 )W^ l \x,p) , (6) 

where = p^ + ^d x ii and the partial derivative with respect to x on the right-hand side acts only 
on the current jO. If it were not for the V-A factor on the right-hand side, the structure of these 
equations would be analogous to Eq. (8) of Ref. |^] and could be analyzed accordingly. 

In order to proceed here, we employ the decomposition of the spinor Wigner functions Ol : 



w (l) = jr(l) + i7 5p(Q + yi V (i) + yy^jf) + 1^5^ , 
i.e., in terms of scalar, pseudoscalar, vector, axial vector, and antisymmetric tensor components: 



A il 



\x,p) = -tr W^(x,p) , 

\x,p) = — -i tr ^W^ix^p) , 

\x,p) = -tr 7 M iy (/) (x,p) , 

\x,p) = -tr 757^W r(i) (a;,p) , 

,)(a;,p) = -tr cr^iy (/) (x,p) , 



which are real functions. Thus, for example, we obtain: 



d 4 p [<®V$\x,p) + <$A®{x t p)) 



(0 



(7) 

(8) 
(9) 
(10) 

(11) 
(12) 

(13) 



using Eqs. (|5|) and (|i~0|), (|Tl|). Only these (axial) vector currents couple the transport equations 
We introduce an abbreviation for the shift operator appearing in Eqs. (ra), 



J, 



(0 



iti 



exp 



• dpWW = I cos (o^ • d p) - ^infe • fit)] JW(x) = 7$ - iZj) , 



(14) 



where d x acts only on the current jji \ as before. Then, making use of the commutation and trace 
relations of the 7-matrices, we decompose Eqs. @ in terms of the Wigner function components, 
Eqs. (p|)-(|l2|). Thus we obtain the set of coupled equations: 



(15) 
(16) 
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J {V)X (Wc{M) - U$V®) ~ i$S% - \cfe» Xw ,S^') , (17) 
J {V)X (g^V® ~ + \$e» Xvv , S W + , (18) 



KK.VP - K V V$>) - + m«S« = J<?» (^ A (c»V« + cJU^), - . . . ^ ] 

-^xA^A^ + c^Y') , (19) 

with as defined after Eq. (^). 

As is well known from other cases [11, [l2|, [l3|], the real and imaginary parts of these coupled 
equations can be separated and eventually will thus lead to the proper transport equations in phase 
space and the generalizations of the mass-shell constraint, cf. Section 2.2 . 

Furthermore, we observe that the left-hand sides of Eqs. (|l~5|)-(|T9|) formally coincide with Eqs. (5.7)- 
(5.11) of Ref. jyj]. There, however, the corresponding operator for electrically charged particles 
necessarily incorporates the effects of the Lorentz force in the external field (Hartree) approximation. 

Due to the linearity of the Dirac equation with respect to the weak and electromagnetic interaction 
terms, i.e. with the derivative in Eq. (|2|) replaced according to the minimal coupling rule, dj 1 — > 



dj 1 + ieA^(x), it is straighforward to incorporate the electromagnetic interaction into Eqs. (15)-(19) for 



the electron-positron case. Making use of the earlier QED results, this is achieved by the substitution: 

K* = IP + ^V^ , (20) 

V m = d£-ej {\d x -d p )F^d p » , (21) 
IP = v ^ - e \^\d x -d v )F^d v , , (22) 



where jo and ji are the conventional spherical Bessel functions, cf. Eqs. (4.19)-(4.21) of Ref. [11]; the 
derivatives d x in their arguments act only on the electromagnetic field strength tensor entering here, 
F^ v (x) = d x u A v {x) — d x u A^(x). Our convention is that e denotes the electron charge. 

With the electromagnetic fields incorporated, we also need to include the Maxwell equation, 

d^(x) = J» m (x) =etrj d 4 p ^W^(x,p) , (23) 

which consistently determines F^" in terms of the electromagnetic four-current Jem- However, an 
important remark is in order here. Together with Eqs. (|20|)-(p2"|) also the definition of the Wigner 
function (B) has to be modified. In order to preserve the gauge covariance of the equations, one has 



to include an appropriate electromagnetic phase factor ('Schwinger string') [11, 12|. Since it will not 
appear explicitly in any of our further derivations or applications, it may presently suffice to keep this 
in mind. 

This completes the derivation of the coupled transport equations for a system of electrons, neutri- 
nos, and electromagnetic fields in accordance with the Standard Model and in the collisionless (Vlasov) 
limit. 
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2.2 The Semiclassical Limit 

Our aim in this section is to extract the relevant semiclassical equations from the quantum transport 
equations which we obtained in the previous section, Eqs. (p^)-(|l9|) in particular. Taking the explicit 
^-dependence into account, which enters through the definitions of the shift and kinetic operators in 
Eqs. (14) and (20)-(22) respectively, it becomes obvious how to expand the equations in powers of 
h. Since the leading terms of the real and imaginary parts of the equations start out with different 



powers, it is useful to separate them, similarly to what was previously done 11, 12, 13 1. 

Furthermore, we presently simplify the set of equations by assuming a spin saturated electron- 
positron plasma, i.e. without the spin polarization effects which may be induced by strong magnetic 
fields, for example. Thus, for the e + e~ plasma, we have no pseudoscalar or axial vector densities, cf. 
Eqs.©-©. 

Also, the Standard Model neutrino-antineutrino system consists strictly only of left-handed neutri- 
nos vl and a right-handed antineutrinos isr, if we appropriately neglect here their tiny (possibly finite) 
masses. In this case, as we show in the Appendix, only the equal vector and axial vector densities 
contribute to the neutrino Wigner function, while all other densities vanish in the massless limit. 

These approximations serve as a working hypothesis for our study of the collective modes and 
their (in)stability in a supernova environment in Section 3 . Eventually, however, the analysis of the 
complete coupled set of equations (|l5|)-([l9|) and (f25|) should be performed, considering the presence 
or generation of strong magnetic fields during supernova explosions or other astrophysical processes 
1 15, 16, |l7]] (and references therein). 



2.2.1 The Semiclassical e + e Transport Equations 

Implementing = = (spin saturation) and separating the real and imaginary parts of 

Eqs. (15)-(19) with the help of Eqs. fli~4| ) and (p0[), we obtain for the e + e~ plasma the set of equations: 

II • V (e) - m (e) ^ (e) = , (24) 

ftV-V (e) = , (25) 

= T& v) -V (e) , (26) 

= xH.y (e) , (27) 

n/ (e) + ^V"5g-m( e )vP = n^\g, x c v ^ -\c A e^S^')-c v l^ x S^ , (28) 

^V^W-ITCW = -l^\g, x c v ^-±c A e, Xuu ,S^)-c v n^ x S { ;l , (29) 

^wIlV^' = -KM x (g^c A ?* e) -~cve^sV^) + c A tf v »S$ , (30) 

je^^sW = jW\ 5M CA^ (e) -^cve^ w <S( e )^') + c^W A 45 , (31) 



\(y^-V v V^)+m^S$ = ^(iMv^-lMv^-^vv^V^', (32) 
- n„V< e > = c^KMyW-^V^ + ce^^Vf^' , (33) 
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where the constraints (|2q ) and ([27]), which resulted from Eq. (|16|) , where taken into account in Eqs. (£ 
and (25), which resulted from Eq. (15); we used Cy A = cy,A , cf. Eqs. (|lj), @. 

We proceed to evaluate the limit ft — > of the above system of equations. To begin with, we obtain 
from Eqs. (JHJ) that ftj? = jjP + 0{fi 2 ) and ijP = \d x ■ d p jf + 0(ft 3 ) and from Eqs. (0), @ that 
IF = p* 1 + <3(ft 2 ) and W = S/ _ ei^o^ + 0(ft 2 ). 

Then, first of all, the vector density can formally be calculated from Eq. (Eq): 



+£(0," - eF^fl^M + c^flS, • <V HA 4 6 a) + °(^) « ( 34 ) 

where the right-hand side is to be evaluated consistently to first order; we recall that d x acts only on 
in the last term. Similarly, we obtain from Eq. (32): 

- M- cAe ^' ji " ,Xv """ 

+\ [(9* - eF„ A 8/)V<=) -...„„„] + cv\a* ■ 9 P (J<"»VW - J<">V< C >)) + 0(R 2 ) , (35) 

with a contribution at O(fi ) in the absence of a pseudo- vector (or -scalar) density, in distinction to 
the QED case of Ref. |11|| . Taking the limit ft —* 0, we solve Eqs. (|34|)-(|3~5"|) in terms of the scalar 
density or, rather, the modified scalar density, 

f(e) (x p) = F {e) (x,p) 



The results are: 

VjS° = ^(Pm-^W b) , (37) 

5 # = ^VA^-^V WA / (e) , (38) 

where we made use of the constraint (p6|), i.e. • V = for ft — > 0, and conveniently added a term 
on the right-hand side of Eq. (|38|) which vanishes identically. Thus, we find that in the semiclassical 
limit the spinor Wigner functions for the spin saturated system are completely determined by the 
scalar density, cf . Eq. (0) . 

Next, using Eqs. (|36[) and (|37|), the Eq. (pBj) yields a transport equation for the scalar density: 



(0/ - eF^dyXpp - <vJM)/W = ((p„ - cvJjrtydf - eF^dy) - c v (d x ■ jM)) /» = , (39) 
i.e. in the limit ft — > . Similarly, we obtain from Eq. ( |24| ) together with Eq. (^) a constraint equation: 
((p - J^) 2 - m( £ ) 2 (1 + (c A /m {e) ) 2 J(") • jM)) /M = , (40) 



where we also used the constraint (26) in the form: 

(p-cyJ^)- jM/( e )=0 , (41) 

which is appropriate in this limit. 
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Clearly, the Eq. (|40| ) demonstrates that it is the kinetic momentum, 

k^^p,-c v J^{x) , (42) 

which should be related to a classical mass-shell constraint. Therefore, we redefine the scalar density 
as a function of the kinetic momentum k, 

f^(x,p)=f^(x,k^ + cvJ M )^f {e) (x,k) , (43) 

instead of the canonical momentum p; we will omit the superscript from , since it is identical to 
the one of in the respective equations. This implies: 

W (6) I P = W (e) l* - cv(dM v) )9kf® ■ (44) 
For the redefined variable and scalar density function, the e + e~ mass-shell constraint follows: 

(V - 2 - cjjM • jW)) /M = , (45) 



instead of Eq. fl4C|). Furthermore, we finally obtain from Eq. ( |39| ) the Vlasov type transport equation 
for the scalar e + e~ density: 

(k-d x - c v k^J^>)d k u - c v (d x • jM) - ek^dv) 

= (fc.^-^M^^-^jM^ + eFH^))/^ = 0, (46) 
rewriting and using here the appropriate leading order in K form of the constraint (j2j 

(^JM^)^(P^ - cy JM)/W = (A; M (^jM^a fe , + (d x ■ jM)) /( e ) = . (47) 
In particular, we also employed Eq. ([37]), rewritten now simply as: 

V« = igy/M . (48) 

We observe that the weak current-current interaction leads to an antisymmetric tensor coupling in 
the transport equation (fHf), which is analogous to the electromagnetic field strength coupling. 

Furthermore, we remark that there are remaining equations of the set (p4|)-([33j) which we did not 
consider here, since the dynamics can be represented completely in terms of the scalar density J-^ e \ 
recall Eqs. (|36|)-(|3"8|). Similarly as in the QED case of Ref. they could be shown to be satisfied 
identically to leading order in the ^-expansion, which we do not pursue here. 

2.2.2 The Semiclassical vlvr Transport Equations and Currents 

For approximately massless neutrinos, with vjt = and T = V = = (see Appendix), we 
obtain a much simpler set of equations from Eqs. (|l5|)-([l9|): 



(K-2J^)- VM = , (49) 
i[(K-2J^)^-...^]-e^ Xl/ ,(K-2J^) x V^>' = 0, (50) 
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using Cy + = 2 ; here = p^ + ^d x » . Separating real and imaginary parts, we expand the 
resulting equations in powers of % : 

(p-2jW)-V M + 0(h 2 ) = , (51) 
(d x + 2d x -d p J {e) ) ■ V M + 0(h 2 ) = , (52) 
[{p-2J^),V^ -...^] + 0(h) = 0, (53) 

e^ Xv ,{p-2J^) x V^' +0(h) = 0, (54) 

thus proceeding similarly as in the case of the electron-positron plasma up to this point. 

However, now it is obvious that the following Ansatz immediately solves Eqs. (53) and ((53): 



vM^p-2j( e V» , (55) 

with a scalar function / of the phase space variables x,p. Furthermore, to leading order in H , it 
converts Eq. (|49| ) into the mass-shell constraint: 

(p-2j( £ )) 2 / H = , (56) 

which demonstrates that it is the kinetic momentum, k^=p^ — 2J^ (x) , which is to be on-shell here. 

Performing analogous steps as in Eqs. d42|)-(|4"4]) before, redefining f( u \x,p) = f( u \x, k) in partic- 
ular, we obtain directly the vl^r mass-shell constraint: 

k 2 f^ = , (57) 

and from Eq. (j52^) the Vlasov type transport equation for the vl&r density function: 

(k-d x - 2k^J^ - d»J^]d k ») /M = , (58) 



which may be compared to Eqs. (fi5| ) and (46) of the electron-positron plasma. 

In order to complete the set of coupled classical transport and constraint equations (Bq), (|46"|), 



(57), and (|58|), we have to reconsider the four-currents entering here and into the Maxwell equation 
(||) in the limit % -> . 

Implementing the spin saturation, in particular A.^ = , and using Eqs. ( J37| ) and (|43|), we obtain 
the weak e + e~ current: 

(*) = / d*p (p„ - c v J^(x))p\x,p) = 4^ / d*k ^f ie \x,k) , (59) 

cf. Eq. (|l^), with Cy = cy . Similarly, the electromagnetic current assumes the form: 

j£ m (x)=4e (d A k^r(x,k) , (60) 
J m e 

cf. Eq. (^). Finally, the weak vl^r current is: 

=8^|/ d 4 p (p„ - 2j(;\x))f^(x,p) = 8^J d 4 fc k,J^(x,k) , (61) 

using once more = A^ , Cy + = 2 , as well as Eq. (|55|) . 
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The closed set of four coupled mass-shell constraint and transport equations, with the currents 
determined by scalar (density) functions, along with the Maxwell equation present the final result 
of our derivation of the semiclassical nonequilibrium transport theory of neutrinos and electrons. It 
incorporates their antiparticles as discussed in more detail for the QED case in Ref. [11], as well as 
electromagnetic fields, assuming an e + e~-spin saturated system in the mean field dominated regime. 

We note that the structure of our final closed set of equations could essentially be anticipated from 
purely classical kinetic theory considerations, as previously observed || |j~2[ . On the other hand, for 
the study of spin-polarization or strong magnetic field effects and higher order quantum corrections, 
we must go back to our previous set of Eqs. (||)-(|l|) of Section 2.1 . 



3 Linear Response Analysis and (Un) Stable Collective Modes 

Presently, we apply the transport theory of Section 2 , in order to derive the semiclassical dispersion 
relations of collective modes of a neutrino-electron system in general (Section 3.1). 

In Section 3.2 we specialize our results for the Type II supernova scenario. The relevant distribution 
functions are introduced in Section 3.2.1 and the necessary response functions calculated in Section 
3.2.2. 

In Section 3.2.3 we evaluate the dispersion relations for various collective modes. We determine 
their (in)stability properties in the neutrino 'beam' plus electron 'plasma sphere' system formed during 
a Type II supernova explosion. The final Section 3.2.4 is devoted to a discussion of the validity of the 
approximations used. 



3.1 The Linear Response Theory for Neutrino-Electron Systems 

The behavior of collective modes, in particular, the onset of instabilities, is determined by the evolution 
of small perturbations of a generic set of stationary distribution functions, which may be caused by 
scattering interactions, for example. Therefore, we write the scalar density distributions in the form: 

f (l \x,k) = fl\k)+5f^(x,k) , (62) 

where fg denotes the assumed homogeneous four-momentum dependent solutions of the mass-shell 
and transport equations and 6f® an initially small perturbation. This assumption of homogeneity 
greatly simplifies the subsequent analysis and describes a sufficiently large 'free-streaming' electron- 
neutrino system. 

The weak currents determined by fjp, cf. Eqs. (|59| ) and (|6"T|), are homogeneous and the 
antisymmetric tensors which enter the transport equations ( |46|) and fl58|) , 

G (i)^ _ c (i)[ d »j(i> _ a^jo>] ; ( 6 3) 

vanish in this case; from here on = cy and c^ 1 = 2 . Furthermore, assuming an isotropic on-shell 
electron-positron distribution, 

f^\k) = 8[k 2 -m^ 2 -cp^ ■ J^]f^(k°,\k\) , (64) 

cf. Eq. ([£^), it follows that the corresponding electromagnetic four-current ( |60|) vanishes, if we addi- 
tionally assume a neutralizing background charge or approximately equal densities of electrons and 
positrons, depending on the circumstances. Consistently we set F^ u = , i.e., considering a situation 
without external electromagnetic fields. 

Indeed, then, the initial on-shell distributions f^\k) and fg (k) are stationary in the absence of 
collisions. They will be further specified shortly. 
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Linearizing the transport equations ( ^6|) and (|58l ) with respect to the small perturbations Sf"\ we 
obtain for the electrons: 



ik • q 5f^(q, k) + ^ (SG^(q) + e6F^(qj) d k ,f^ e) (k) = , 



(65) 



where we introduced the Fourier transform for any function g of the space-time coordinates, g(x) = 
(27r) — 4 / d 4 gexp(— iq-x)g{q). Here 5G^ and 5F denote the weak and electromagnetic tensors induced 
by the perturbations Sf( v ' and 6f( e \ respectively. From Eq. (63) we obtain: 

5G^ v {q) = -ic il) {q»g uX - q u g pX )8J { p . (66) 

Furthermore, solving the Maxwell equation ( |23| ) with the retarded boundary condition (damping in 
the infinite past), we obtain: 



eSF^iq) 



ic 



q* + leq 



o(«"«)^< 



em 



V2e 2 



1 



cvGf q + ieq' 



(67) 



where e — ► + , and where we used Eqs. (59) and (|60|), in order to express the (conserved) electromag- 
netic current fluctuation in terms of its weak counterpart. 

Implementing the retarded boundary condition, i.e. the 'Landau prescription' ||], the electron 
transport equation (|65|) is solved by: 



8f( e \q,k) 



k^ (dG^(q)+eSF^(q))d k , 



f { s\k) 



—i(k ■ q + iek°) 

Similarly, the perturbation of the stationary neutrino distribution is determined by: 

—i(k • q + iek u ) 



(68) 



(69) 



Obviously, Eqs. ( pq ) and (6£) are coupled to each other via Eqs. (66). 
We proceed by introducing the response functions: 



M^ xp (q) 



d 4 k 



1 



(k-qd k p -k»q-d k )f$\k) 



( l=e ) k-q + iek° 



m 



(70) 



which will be calculated for specific choices of the stationary distributions fjj' shortly; the factor 
l/mS l=e ^ is meant to apply only in the e + e~ case and to be replaced by 1 for the (approximately) 
massless vlvr case. 

Making use of the response functions, we multiply Eqs. (|6g|) and (6£) by the appropriate factors, 
cf. Eqs. ( |59|) and (]6l|), and integrate over d 4 fc, in order to obtain a closed set of algebraic equations: 



8jV\q) = M^(q)(^f5J^(q) 



q 2 + ieq' 



(71) 



5J^ X (q) 



4G F 
V2 



M^ X P{q)5J ( f\q) , 



(72) 



where the tensor fluctuations 5G^ and 5F were eliminated with the help of Eqs. ( |66| ) and (|67|), 
respectively. Inserting the second into the first equation, the final result is: 



M^(q)dJ^(q) = 



7 x P + M^ Xa (q) 



q 2 + ieq 



gp - 2c 2 GharM^(q) 



8A e \q)=0. (73) 
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The solvability condition of this vector equation determines the dispersion relation for the perturba- 
tions of the stationary electron-positron distribution: 

Det M {e) {q) = , (74) 

where Ai^ is a 4 x 4-matrix in the Lorentz indices. Analogously one obtains the dispersion relation 
for the neutrino case, which we do not pursue. 

A final remark is in order here. From the structure of Eq. (|73|), particularly the generically small 
weak coupling term compared with the electromagnetic one, it is natural to expect that the neutrinos 
can only influence the resulting dispersion relations noticeably, if their response function shows rather 
singular behavior. Furthermore, in this case, weak and electromagnetic interactions presumably will 
mix in the corresponding collective modes, due to the products involved, for example, in Eqs. ( |7"3| ) and 
(J74|). This will be studied in the following sections with the application to a supernova scenario. 



3.2 The Supernova Two-Stream Scenario 

The above results are fairly general and need to be specialized according to the physical nature of 
the stationary distributions as well as of their potentially unstable perturbations. We shall now study 
the idealized situation where an electrically neutral finite temperature electron-positron plasma is 
hit by a neutrino-antineutrino beam. (Anti-)neutrinos are radiated from the neutrino sphere, move 
approximately radially outwards, and interact collectively with the electron-positron plasma sphere 
forming the 'radiation bubble' in Bethe's supernova scenario 0. As before, we derive our results in 
the collisionless limit. 

Typically, it is assumed that a short-lived, but intense neutrino flux (3 x 10 29 W/cm , total 
integrated luminosity up to the order of several 10 53 erg ) with an approximately thermal spectrum 
corresponding to a temperature T v ~ 1 . . . 10 MeV is released from the collapsing core and interacts at 
a distance of about 30 . . . 300 km from the center with the surrounding moderately relativistic electron 
plasma of (charge) density n e ~ 10 30 cm -3 and temperature T e ~ 0.5 MeV; here the uncertainties 
mostly reflect differing scenarios considered in this context j|, ||, |||, [j| [|. We will study a corresponding 
set of parameters, following the discussion of the radiation bubble by Bethe ||. 

For the above (optimistic charge) density and temperature, the electron-positron plasma is non- 
degenerate, with an estimated chemical potential fi e < ttT . This is indirectly supported by Bethe's 
results, see in particular Sections VI. E-G of Ref. [||, which demonstrate the dilute character of matter 
in the radiation bubble - the energy or entropy density of the 'radiation' (i.e. of photons plus pairwise 
produced electrons and positrons) is more than a factor 10 2 higher than that of nucleons in the bubble. 
Therefore, any background charge contamination by protons must be small here, and correspondingly 
the net charge of electrons over positrons neutralizing the plasma. Hence we may neglect the finite 
electron chemical potential in a first approximation. 



3.2.1 Distribution Functions 

The following stationary electron-positron distribution will now be considered, cf. Eq. (]& 

f( e \k) = (2Tr)- 3 m^5[k 2 - 2 - cjj^ • jM] (e(k°)F(k° /T e ) + e(-k°)F{-k° /T e )\ , (75) 

where F(x) = (e x + 1) _1 , and where T e denotes their temperature. When jM = , Eq. (75) describes 
the e + e~ blackbody radiation (omitting the vacuum contribution). We remark that antiparticles are 
represented as fermions with negative four-momentum here fl~T|| , 

Concerning the emission from the neutrino sphere, we neglect its collective flow relative to the 
electron-positron plasma sphere, or vice versa. However, it is important to incorporate the dilution 
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and angular squeezing effects due to the spherical geometry. Thus, we assume the following stationary 
(approximately massless) neutrino-antineutrino distribution: 

f£\k) = i(2vr)- 3 5[fc 2 ] (Q(k°)e(9 max -6 HJ :)F([k - ^]/T u ) + ... k,^ - -k-^) , (76) 

where [i v denotes their chemical potential and only one or vr) spin state is taken into account. 
The additional ©-function, implementing the radial ('outward') unit vector ft, accounts for the finite 
opening angle 6^ r between neutrino momenta and the radial direction. The maximal opening angle 
is determined by sin Oyyiax — Rjv, where r denotes the distance from the center of the neutrino sphere 
of radius R. The usual dilution factor, d = (R/r) 2 = sin 2 9 max , does not appear explicitly, but is 
recovered in the calculation of, for example, the energy flux from the neutrino sphere based on f^. 

We remark that the neutrino distribution is not necessarily uniform within the cone defined by 
9 max- It may vary considerably, depending on the emission characteristics of the neutrino sphere. 
Thus, the distribution of Eq. (^) may represent an opening angle average; the (ir)relevance of the 
sharp 0-function cut-off will be discussed in the final subsection 3.2.4 . Furthermore, it is not a global 
solution of the spherical free-streaming problem. However, our approximate spatially homogeneous 
distribution, with the parametric dependence on R/r, is sufficient for the study of collective modes 
with a characteristic wavelength very much less than R, even though we are interested in the long- 
wavelength limit with respect to the microscopic scales. 

We omit the n.- and r-neutrinos at present which have a considerably weaker effective coupling, 
see Eq. (|73|) together with the remarks after Eq. (|l|); using sin 2 6\y ~ 0.23, we have cy « 0.96 (—0.04) 
for the electron (fi, r-) neutrinos. Furthermore, their chemical potential vanishes. For the electron 
(anti-)neutrinos, which carry about 4/10 of the total neutrino energy flux, we adopt Bethe's estimate 
which yields r\ v = [i v jT v = 0.29 p. 

Next, we proceed to calculate the energy-momentum tensor, similarly as in Ref. [11], for a stationary 
free neutrino-antineutrino distribution: 

T$ (x) = tr J d A k k vlll W [ ^ (x, k)=Af d 4 k fc^vM (x, k) , (77) 

i.e. in terms of the vector density, cf. Eq. (|To|). Employing Eq. ( |55|) and projecting on the 'outward' 
momentum direction, we obtain the electron- z^i>r energy flux corresponding to the homogeneous 
equilibrium distribution (|76|): 



T (v)0i n i = 4 / d A k k o k z f M (k) = JL / dkk 3 _ ^yrpj + F{[k + ^j/r,)) 







8vr 2 

= ^/fl + 4^ + 4^=1-^ , (78) 
480 " V 7vr 2 lu 7vr 4 ,u ) 10 irR 2 ' V ; 

choosing n = (0,0,1), and where L denotes the total neutrino-plus-antineutrino luminosity. The 
integral is evaluated exactly with the help of a formula from Ref. [Jig]. Correction terms involving 
powers of m^/(r] u T u ) would be completely negligible for temperatures in the MeV range and a 
typical neutrino mass (much) less than 1 eV. The last equality in Eq. ( fTq) provides the relation between 
temperature and radius of the neutrino sphere, given its luminosity L 1|. 



Similarly, we obtain from Eq. (61) by direct calculation: 

t3 



jW, = 8 ^= J d*k k»ff\k) = S|f§. + JLtf) e , (79) 

where ^ = (2/[l + cos max ] ; 0, 0, 1) . As expected, the neutrino current components are very small, 
since for temperatures of about 10 MeV we have GfT u 3 ~ 10 -8 MeV. Therefore, the corresponding 
term oc 

jM . jM i n the expression for the stationary electron-positron distribution, cf. Eq. ([75]), can 
be safely neglected henceforth. 
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3.2.2 Response Functions 

After specifying the unperturbed stationary electron and neutrino distributions, and fg respec- 
tively, we calculate the response functions M^ xp and M^ xp defined in Eq. (f70|). For the following 
calculations it is convenient to perform a partial integration, which yields: 



M^ xp (q) 



d 4 /t 



JS 



(k) 



in 



(l=e) 



_ g x P+ <l XkP + <l Pk 



A 



k x k p 



(80) 



k-q + iek° (k-q + iek°) 2 J ' 

0. Thus, the current fluctuations are 



which is now obviously symmetric and transverse, q\M^ Xp (q) 
properly conserved, q\SJ lx (q) = 0, cf. Eqs. (ffl])-(|7^). 

Beginning with the electron case, the calculation is facilitated by recalling that the distribution 
function fj?\ Eq. (|75|), is isotropic with respect to the three- momentum components. Therefore, the 
spatial part of the response function can be decomposed into a transverse and a longitudinal part, 



M^ ij (q) = [5 ij 



n2 



M, 



'{q) + —*-Ml>(q) 



Defining the electric (Debye) screening mass, 

4p2 poo i 

m 2 D = — dk kF(k/T e ) = -e 2 T e 2 , 

TT Z Jo O 

the results of a standard calculation for the electron-positron response function are: 



2 m( £ )°V, 

M^ 0i (q°, 



m 2 D 



r i« 




q + oj 




L 2q 




q — to 



iirO(q — to] 



M^ i0 (q° 



. M (e)00 (g )( 



(81) 

(82) 

(83) 
(84) 



M ( L e \q°,q) = ^M^°(q ,q) 



e 2 M<t e \q°,q) 



;m D - 



uJ- 





Mn 


q + oj 




. 2\q ujJ 




q — ijj 



inQ(q 



•>) 



(85) 



(86) 



where we implemented the relativistic limit and neglected correction terms in powers of m e /T e ; here 
we simplified the notation by introducing uj = and q = \q\. We note the appearance of the 
imaginary parts which, in general, are responsible for Landau damping |], [ll|. These kinetic theory 
results completely agree with the perturbative one-loop evaluation of the QED polarization tensor in 



the high-temperature limit [20], a correspondence which was has been observed in many other cases, 
e.g., 0. 

In a fully realistic calculation, corrections due to the finite ratio of m e /T e ~ 1 should also be 
considered. We neglect them in our present work, since the appropriate electron temperature is 
not precisely known in this context. Furthermore, unfortunately, this would necessitate numerical 
calculations where the transparency of the analytical results presented here would be lost. On the one 
hand, it seems unlikely that the additional mass scale can qualitatively change any of our conclusions, 
since it is well separated from all the plasma scales entering in the following. However, for particular 
effects, e.g. proper Landau damping || 19], a finite mass may be crucial (cf. footnote following the 
discussion after Eq. (112)). Thus, proper 'neutrino Landau damping' has been discussed in more detail 
by Silva et al. recently ||. 

Next, we turn to the calculation of the neutrino response function M^ x f '. It is more involved due 
to the preferred direction of propagation, which enters here through the dependence of the stationary 
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distribution (|7q ) on the 'radial' unit vector n. In order to facilitate our task, we consider two cases 
separately, depending on the orientation of the wave vector q with respect to n: q \\ n (Case I) and 
q _L n (Case II). We recall that q determines the direction of propagation of the collective excitations 
of the electron-positron plasma, especially in the presence of the neutrino flux. 

Case I. Here we expect a response function with a formal structure generalizing the familiar results 
of Eqs. (p^)-(p6|), since the geometry determining the essential angular integrations is identical to the 
previous case. Therefore, a tensor decomposition into transverse and longitudinal parts analogous to 
Eq. (^) still applies. However, the maximal opening angle 9 ma x between neutrino momenta and the 
'radial' direction limits the azimuthal angle 9, e.g. in the vector decomposition, 



k = qq 1 k cos 9 + kj 



(87) 



with k = \k\, which is conveniently employed after converting the integral of Eq. ( |80| ) to the corre- 
sponding threedimensional (on-shell) form. 

Furthermore, instead of the Debye mass of Eq. ([32"1), we introduce the weak thermal mass: 



m,. 



2c 

7T 



poo i / 3 

J dk k (F([k - fx v ]/T v ) + F([k + Hu]/T v )) = -4T U 2 \1 + -^j 



which takes the finite chemical potential of the neutrinos into account. Here we applied the ultrarel- 
ativistic limit discussed before, as well as the appropriate integral formula from Ref. fl8|] . 
Then, we obtain the components of the neutrino- antineutrino response function (q\\ n): 



2c 2 M^ 00 (q°,q) 



M^ 0i (q°,q) 
Mi U \q°,q) 
2c 2 M^\q\q) 




iir@(q - q°) + iirQ(zq - q° 



M^ 00 (q°,q) 



(89) 

(90) 
(91) 



Z ( l + zq(l-q 2 /uj 2 ) 




zq 
zq - 



iTT&(q - q°) + iTT@(zq - q°) 



,(92) 



where z = cos9 rnax and u> 2 = (q ) 2 - Indeed, for z — > —1, i.e. without restriction on the opening angle, 
we recover the formal structure of Eqs. (|33|)-(|8"6|), while for z — > 1 the response function vanishes. 

Furthermore, we observe that for a finite opening angle (1 > z > — 1) the Landau damping 
imaginary parts are limited to the region q > q° > zq and that at the resonance frequency q° = On = zq 
the response function has additional singularities, which are absent for z = —1. In the present case, 
with q || n, the longitudinal as well as the transverse components, and respectively, are 

affected. 

Case II. In this case, with q _L n, we introduce a third unit vector e, perpendicular to the other 
two vectors, in order to decompose the momentum vector for the threedimensional response function 
integral, 



k = k(n cos 9 + qq sin 9 cos 4> + e sin 9 sin < 



(93) 
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with the azimuthal and polar angles 9 and (ft, respectively, such that k-q = kq sin 9 cos (ft. The resulting 
angular integrations can all be done analytically in the appropriate ultrarelativistic limit, either by 
elementary or contour integration techniques. 

Due to the symmetry properties of the required integrals, presently it turns out to be useful to 
decompose the spatial part of the response function as follows: 



M^ ij (q) = [5 ij 



VbTI? 



M$\q) + n % n?M%>(q) + ^T M lJ(q) + 



Ay) 



(")/ 



n l q J + n J q l (iy) 



M^'(q) . (94) 



All other terms which could arise vanish identically, since the corresponding polar angle integration 
comprises an odd function. 

In order to check the ensuing lengthy calculations, we also evaluated independently the integrals 
resulting from the transversality condition mentioned after Eq. (|80|), q\M^ Xp (q) = 0, as well as from 

the trace M^ X (q). These results we compared with what is obtained using the calculated components 
of the response function in the following. In fact, this procedure leads to considerable simplifications. 

Then, for < 9 max < tt/2, i.e. < z < 1, we finally obtain these components of the neutrino- 
antineutrino response function (gin): 



2c^M^ 00 (q°,q) 
M^ 0i (q°,q) 
24MP{q\q) 

24MtKq\q) 
M%(q°,q) 



vn„. 



lib 1 



z z 



o 



+ 



= =F-— ln(- ^ 

tz 2 1 zq+ ^Jqi 



M^\q\q) = ^Mt\q ,q) + ^M^(q ,q) 



1 2 (^ 2 

1 n 



1 



1 



Tin. 



l-z±( 



q + <t 



2 2 \ q 

(qW 



q°) zq + ^/q* 



1 1 



q ° ln(- q + Q ° 



q 



' zq + 



^M^ 00 (q°,q) 



,1 2 9° 



0\2 



-(1-(9 U ) 



Q\-l 



0\2 



l)(l-5°/V5) 



(95) 
(96) 
(97) 

(98) 
(99) 
(100) 

(101) 



where we introduced the abbreviation, 

^(g°) 2 -? 2 (l-^ 2 ) , 
with q = \q\ and z = cosO max , as before. Several qualifying remarks are in order here: 

• For later convenience we did not separate real and imaginary parts in Eqs. (p5|)-( p!00| ) which are 
valid for complex q° = uj + ij, provided the imaginary part here is sufficiently small, I7I <C \lo\, 
or infinitesimal. 

• Either the upper or the lower signs have to chosen consistently in Eqs. (|95|)-( [TO0| ) according to 
the following rules (0 < 9 max < 7r /2): 



7 > and uj > 
7 < and uj < 



upper signs ; 7 > and uj < —qsm9 max 
lower signs ; 7 < and uj > q sin 9 max = 



> lower signs ; 
upper signs . (102) 



They are due to the (angular) contour integrations, which result in different contributions ac- 
cording to the listed rules. 
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We do not report the results for u in the intervals which are excluded in ( |102| ), since the azimuthal 
angle integrations have to be split in this case, yielding even more complicated expressions. 

Obviously, the response function has additional square-root singularities at the resonance frequen- 
cies iO = = ±q sin Qrnaxi *y — ► 0. The transverse component Af^p \ however, is not affected in the 
present Case II (g 1 n). 

This completes the calculation of the response functions for the model distributions discussed in 
the previous subsection. 



3.2.3 Dispersion Relations, Collective Modes and Instabilities 

It is useful to begin the study of the dispersion relations following from Eqs. (|73|)-(|7^) with the case 
of the electromagnetically interacting electron-positron plasma, i.e. with the weak interaction term 



tx Gp in Eq. (73) switched off. 

Considering separately transverse ('T") and longitudinal ('17) current fluctuations, i.e. 5J^ e \q) _L q 
and 5j( e \q) \\ q, respectively, Eq. ( ]74| ) yields two equations determining the corresponding dispersion 
relations: 

T : (l-e 2 M^ 1 /(q 2 +ieq°)Y = , (103) 



T 

L: 1 - e 2 M^ 00 /q 2 = , (104) 

with the plasma response functions of Eqs. (|33"|)-(|36|), and where the decomposition (81) of the spatial 
part of the response function is especially taken into account. 



The real solutions with u = \q \ > q = \q\ of Eqs. ( 103 ) and ( |104| ), respectively, determine the 
well-known collective transverse and longitudinal plasmon modes |9|, |20(| . In the long- wavelength limit 
(lo ^> q), for example, the explicit solutions are: 

T : co 2 (q) = uj 2 + p 2 , L : ^l(q) = co 2 + p 2 , (105) 

with the plasma frequency loq = \m 2 D , cf. Eq. (^), which characterizes an ultrarelativistic neutral 



3 

plasma. Again this is in agreement with the one-loop calculations of finite temperature field theory 
(T 3> m e ). We remark that beyond the present collisionless approximation these modes naturally 



aquire a finite width [2C 



We now turn to the case of a fully interacting neutrino-antineutrino beam impinging on an electron- 
positron plasma. We remind ourselves of the two limiting cases introduced in the preceding section 
concerning the orientation of the wave vector q with respect to the outward normal vector n, i.e. 
q || n (Case I) and q _L ft (Case II). In both cases, we concentrate on the interesting possibility 
that the weak interaction term might become comparable to the purely electromagnetic term oc e 2 in 
Eq. (f74|). Due to the intrinsic smallness of the weak coupling constant this may happen only, when the 
neutrino-antineutrino response functions become large, close to the singularities found in Eqs. (|8^)- 



(92) or in Eqs. (|95[) - ( |100D . Otherwise the neutrino effects can be treated as small perturbations of 
previous electron-positron plasma results, as we shall see. 

Case I. We observe here that the neutrino-antineutrino response function obtained in Eqs. (|89|)- 
J92] ) has the same tensor structure as the electron-positron one of Eqs. (p3|)-(p6|). Considering the 
product of the two appearing in Eq. @, C a/3 = M {e ^M^ ] 13 , we obtain: 



C oo = (l _ ") M ^ 00 M^ 00 , (106) 
C 01 = C 10 = ^C 00 , (107) 
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q 2 q 2 



C*i = - ( ^ - ] M^M^ + ^C 00 , (108) 



with w = \q°\ and g = Therefore, we may still distinguish transverse (T) and longitudinal (L) 
current fluctuations which do not mix, similarly to the case of a purely electromagnetic plasma. 

In analogy to Eqs. ( |103D and ( |104j ), we thus obtain from Eq. (^|) two equations which now determine 
the neutrino 'beam' electron-positron plasma dispersion relations (q\\ n): 



e 



2 



T: 11- i ^2_ q2 + 24GM^j = , (109) 

L: 1 - [4 + 2c^G |(1 - i^)2 M H 00 ]M( e ) 00 = , (110) 
q 2 q 2 

with q = \q\. 

In the long- wavelength limit (uj 3> g) and to lowest order in Gp we obtain, for example, from 
Eq. ( |109D the equation: 

= u: 2 -q 2 -u;l{l + \{l? + ...) 

A(l-z)+z(l-z 2 ) (l-z 2 ) 2 q 4(1 - z) - 3z(l - z 2 ) 2 q 2 
-u, a ( — - + (-) +...), (Ill) 

with uj = | q° |, z = cos6 max , and where we indicated the neglected higher order terms in q/uj. We 
introduced the dimensionless constant: 

« 2 = T2 G F™1 m l , (H2) 
2e z 



which governs the strength of the neutrino effects. The solution of Eq. ( Ill ) describes the transverse 
plasmon in a neutrino 'beam' electron-positron plasma. 

However, as we anticipated, the smallness of the weak coupling constant makes the influence of 
the neutrino terms completely negligible here. Considering Type II supernova conditions and setting 
T e ~ 1 MeV and T u ~ lOMeV, we find that a 2 ~ 10~ 22 . Omitting the neutrino contribution and 
solving reproduces the first of Eqs. ( |105| ). 

A similar analysis, i.e. for uj > q, applies to Eq. ( |110| ) which the longitudinal plasmon in a neu- 
trino 'beam' electron-positron plasma. Again the neutrinos have a negligible effect under supernova 
conditions. 

We now consider the dispersion relations implicit in Eqs. ( |109| ) and ( |110| ) close to the resonance 
frequency, i.e. uj ~ Oil = zq, which lies in the electron-positron Landau damping regime with < uj < 
q, considering < z < 1 from now on (0 < 6 max < vr/2).0 

In this case, we expect the frequency q , and correspondingly uj, to aquire a finite imaginary part, 
instead of the infinitesimal ie representing the retarded boundary condition ||, cf. Eqs. (|68|), (|69|) 
or (|80|). Therefore, replacing q° + ie — > uj + ij, the 'Landau logarithms' and imaginary parts of 
the calculated response functions have to be reconsidered. We rewrite uj = zq + £, anticipating that 
£ <^ qz, and will use: 

uj + q + i-f 1 + Z . i /-* , i + il v W1 i + il , niQ x 

In — =ln- zvrSign7 + ln(l + ———)- ln(l - — -) , (113) 

uj — q + vy 1 — z q(l + z) q{l — z) 



For the case of an ultrarelativisti c pu re electron-positron plasma in equilibrium it can be shown that no solution, for 



example, of the dispersion equation (104) exists with < ui < q 
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where Sign7 = 7/I7I . This is most appropriate for small £ and 7, reproducing the usual result for 
7 0+. 

Specifically, we reconsider Eq. (110) and take only the dominant singular term oc (zq — q )^ 1 in 
M^) 00 into account, cf. Eq. ©. Thus we obtain more explicitly: 



1 + 



q 2 2 q 2 zq — u> — 27 



1 w + 17 , u) + q + i-j 
1 in 



2 q u — q + «7 



, (114) 



with u; = + £. We recall that a 2 <C 1. 

It is easy to see that for a solution with u> ~ zq the term oc a 2 has to behave qualitatively such 
that (at least) {^,, r f)/q ~ a 2 q 2 /m 2 ~ ) <C 1, particularly in the long-wavelength limit with q 2 /m\) <C 1. 



Consequently, using Eq. ( |H3|) , we expand Eq. ( |114j) up to second order in £/q or j/q. Separating real 
and imaginary parts, it is straightforward to solve the resulting equations. We obtain: 

£ 1 2, , Wl 2 ^z/2) 2 +f 2 (z)+f{z){q/m D ) 2 

- = -a (q/mn) (1 - z ) ^ - ^ + ^— ^ , (115) 

I7I 7T 2/ 4 (1 — Z 2 ) 3 

T = 2° (vrz/2) 2 + [/(*) + ( q /m D ) 2 ] 2 ' (U6) 

neglecting higher order in a 2 corrections and defining: 

/(z) = l-izlni±| . (117) 

These results are consistent with the applied expansions, noting that £ oc Gp/e 2 and 7 oc Gp/e 4 , 
particularly in the long- wavelength limit. 

Recalling u = zq + £, we thus obtain a pair of longitudinal pharon modes ('Type I', i.e. for q\\ n), 
with the real part of the dispersion relation in the long-wavelength limit given by: 

tu(q) = zq+^G$m 2 w (l-z 2 ) 3 q 3 , (118) 

one, a growing and the other, a decaying mode, depending on the sign of 7. 

Analogously, we analyze the transverse dispersion relation to be calculated from Eq. (109) for 
u < q. In this case, we find a pair of transverse pharon modes (Type I) with: 

1 -L 2 (n/m Wl ^3 W 2 ) 2 (l ~ z 2 ) + (z 2 / (1 - z 2 ))g 2 (z) + 2g(z)(q/m D ) 2 

q ~ 4 aW/mDMi Z) {vz/2) 2 + [{z 2 /(l-z 2 ))g(z)+2{q/m D ) 2 ] 2 ' {LU) 



h ' V(g/m p ) S ^ ; F7 pTTT (1 ^w\ I 57 / ^ > ( 120 ) 



? 4 - v-/ (7rz/2 )2 + [ (z 2/ (1 _ Z 2)) 5(z) + 2 (q/m D ) 2 ^ 

and where: 



= 1- Uz- -)ln^f . (121) 
yv ; 2 V z 1-z y ' 

The corresponding real part of the dispersion relation in the long-wavelength limit is: 
M X ^ 22 n ^ 3 (V2) 2 (l-z 2 ) + (z 2 /(l-z 2 ))g 2 (z) 

with an interesting negative sign in front of the second term. We observe that the transverse pharons 
are quite sensitive to the geometry parameter z. 

In particular, in the limit z — » 0, corresponding to a maximally fanned-out 'beam' with 9 max — > ir/2, 
the 'damping constant' 7(g) diverges. In this limit the expansions leading to Eqs. (|119|) and ( |120D 
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clearly break down. This can be studied in more detail, starting again with Eq. (|109| ) and implementing 
q° = £ + i"f, with £,7 <C «?• However, it leads to a nonpropagating mode with frequency of the same 
small order of magnitude as the damping constant, which is physically irrelevant to our study. 

However, under Type II supernova conditions, with a 2 ~ 10 -22 , and recalling that we have 1 — z 2 = 
sin 2 9 m ax = {R/ r ) 2 , in terms of the radius R of the neutrino sphere and the distance r of the electron- 
positron plasma from its center, a typical value may be (R/r) 2 » 0.5 for R ss 30 km [|. Then, for 
q < mo, the above calculations are accurate and we may roughly estimate, for example, the transverse 
pharon damping constant, 7 ~ 10~ 2 a 2 (q/mrj) q- For a pharon wavelength corresponding to q ~ m.o/2 
and an electron temperature T e ~ 1 MeV, this yields a growth/decay length (one e- folding) on the 
order of 10 9 — 10 11 km. A one-percent increase of the collective mode amplitude squared, i.e. of its 
energy, means it would have to run through more than 10 6 km of plasma, which is simply not there. 
The longitudinal mode behaves similarly. 

Clearly, the above estimates are crude and could be improved by folding the results with the 
appropriate distributions, depending on the distance from the supernova core (and time). However, in 
view of the intrinsic weakness of the instabilities, we conclude that it is unlikely that long-wavelength 
Type I (q \\ n) pharon modes play an important role in the outward energy transport processes in 
Type II supernovae. 

Considering the strong momentum dependence of the calculated damping constants, Eqs. ( |1 16[ ) 
and (|120| ), however, the question is raised, whether, at shorter wavelengths, corresponding collective 
modes could become important instead. As we will discuss in the following subsection in more detail, 
in this limit, the presently employed semiclassical transport theory breaks down, necessitating further 
study. 

Case II. We recall that here we have q _L n and proceed as before. However, the product of the 
two response matrices appearing in Eq. @, C afi = M {e ^M ( " ] /3 , has to be recalculated. Taking the 
different tensor structure of M^, according to Eqs. (p3)-(100), into account, we obtain: 



C 



00 



(l_^)M( e )»MM 00 , (123) 



c oi = £9 c oo + n i,9 _5L )M (e)oo M M (124) 
q l q q 

C i0 = «V C 00 _ ^ M (e) M M ; (125) 



« = ^^C 00 + ^(1 - EmW»m 3 m 



C v.. 

q 2 q 2 q q 2 

- (V - nV - ^r)M { T ] + nV'4 1 + , (126) 

with q=\q\. We observe that C"" / C Ufl . 

In the present case, we consider again two different kinds of current fluctuations, when evaluating 
Eq. ([[!): 5j( e \q) lg,n ('Out'), i.e. fluctuations which are perpendicular to the plane spanned by q 
and n, and fluctuations with 5J^ e \q) in this plane ^In 1 ). Thus we obtain the following two equations 
which determine the neutrino 'beam' electron-positron plasma dispersion relations (q _L n): 

° ut : ( X " \ q oy_ q 2 + ^c 2 G 2 M^]M^ = , (127) 
In: (l - 24G$M^M^) + 2c^G|(1 - ^) 2 M^ m ]M^A 
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+(2c^Gj) 2 (2 



(q 



0\2 



(g 



0\2 



-)(M 



M (e) M ( e) 00 = 



(128) 



with q° = uj + £7. We observe that Eg. (|127| ) has the same formal structure as Eq. (109) before. 
Furthermore, if we set Gp to zero, these equations reproduce the transverse and longitudinal electron- 
positron plasmon dispersion equations 



10D and dlQ4[) . 

Guided by our analysis of the transverse plasmon dispersion relation under neutrino flux for q \\ n, 
Eqs. (|109| )-( |112j ), we expect only a negligible perturbative influence of the neutrino interactions in 
Eq. ([L27D, since they are again suppressed by the factor a 2 « 1CT 22 . In particular, the present structure 



of Mf J is a smoothly deformed version of compare Eqs. (plf) and (|9"7|), with no additional 



singularity. Thus, the corresponding collective mode shows no particularly interesting behavior and 
presents a second kind of perturbatively deformed transverse plasmon. 

Finally, we consider Eq. ( |128| ), describing a geometry which resembles the one where two-stream 
instabilities arise in other plasmas ||, ||, 10]. We attempt to find pharon type solutions in the present 



case as well. For this purpose, we take into account the leading root-singular terms, which contribute 
here from Eqs. (pB|), (|98|), and (|100j ). The singularities occur at the resonance frequencies uj = 17^ = 
±gs, as q 2 z = (q ) 2 — q 2 (l — z 2 ) — > 0, where q = \q\ and z = cos 9 max . 

We concentrate on the positive frequency solutions of the dispersion equation ( |128| ) and are partic- 
ularly interested in those with a positive imaginary part, which grow exponentially in time. Therefore, 
we consider q° = u + ry, with uj qs > 0, defining s = sin Offiax- TliGn, Eq. ( |128[ ) ctssumes a slightly 
simpler form: 



= l + m 2 M^ 00 + a 



2 z q° r (?°) 



0\2 



fz q 2 



1] _ l]M(e)00 + lM(f)[l +m ^(e)00]j +0(fl4) {ug) 



with the dimensionless effective coupling constant a 2 , Eq. ( |1 12| ) . The terms of 0(a 4 ) will be neglected 
in the following, since their singularities cancel. Furthermore, we conveniently define: 



ml = m 2 D /q 2 , M^ 00 



l M^ 00 /m 2 D 



M, 



2e 2 MP/m 2 D , 



(130) 



cf. Eqs. (^) and (]8q). Recalling the smallness of a 2 , it is obvious that any interesting solution must 
arise close the resonance frequency Qj_ = qs (0 < s < 1). 

Setting to = qs + £ and assuming £, 7 <C gs, it is useful to expand Eq. (|129|) in terms of the small 



complex quantity k = (£ + i r y)/(2qs). Here we make use of Eq. (113) once more, in order to expand 
the Landau logarithms. Then, expanding to leading order in k, it is straightforward to arrive at the 
'formal solution': 



£ + 27 = qa 4 



s(l — s 



2\3 



128[(l+m2/( s )) 2 + (vr S /2) 2 ] 2 



(/ l ( S )[l+m 2 (/( S )-(m S /2)Sign 7 )]) 2 , (131) 



where: 



h(s) = 2s 2 g{s) + iirs(s 2 - l)Sign 7 + 4(s 2 - l)(/(s) + (ivrs/2)Sign7) 
+m 2 Jf(s) + (ivrs/2)Sign7)(2s 2 5 (s) + ms(s 2 - l)Sign 7 ) . 



(132) 



The functions / and g were defined in Eqs. ( |117j ) and ( |121| ). The appearance of Sign7 on the right-hand 
side restricts the possibility of an explicit solution. 

After some algebra, one obtains a criterion for a solution to exist in the relevant regime (0 < s < 1): 



(_1 + s 2[i + m 2 g g(s)/2]) (vr 2 ™ 2 * 2 + 4/(s)[l + m 2 /(s)]) + 2s 2 g(s)[l + m 2 q f(s)) < , (133) 
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which in the long-wavelength limit (m 2 > 1) can only be fullfilled for sufficiently small s, i.e. sufficiently 
small opening angle of the neutrino momentum distribution. Clearly, taking only the leading terms 
in this limit into account, no solution exists. On the other hand, for mn/q = 2, for example, the 
solvability criterion requires s < 0.426, corresponding to 9 max ps 25°. 

It is obvious, however, from Eq. ( |13l| ) that any pharon ('Type II', i.e. for gin) solution here 
will have £,7 oc a 4 oc Gp, with no particular factors especially enhancing the damping constant. We 
refrain from giving explicitly the not very illuminating lengthy expressions. 

Instead, we conclude that in the supernova environment the growth rate of the presently studied 
Type II pharons is suppressed by an extra factor of a 2 w 10~ 22 , as compared to Type I. Consequently, 
these modes do not contribute at all in this case. 



3.2.4 Discussion 

The detailed calculations in the previous subsections are based on the stationary electron and neutrino 



distribution functions, fg (k) and fg\k) of Eqs. (|75|) and (76), respectively. While f s e \k) is adequate 



for the supernova scenario discussed here, the question arises as to whether the sharp azimuthal angle 
cut-off, present in fg(k), may not cause spurious effects or invalidate our semiclassical transport 
approach. 

In fact, the semiclassical approximation of Section 2.2 is based on the expansion of the full quantum 
transport equations in powers of h, appearing especially in the dimensionless combination Tid x ■ d p , 
cf. Eqs. (H) and (|20|)-(p2|) in Section 2.1. Therefore, a cut-off on a spacelike momentum coordinate, 
corresponding to the angle 9 max between three-momentum and outward normal direction, may produce 
large higher order corrections. These are controlled, however, in the long-wavelength limit. We recall 
that in the derivation of the linear response theory in Section 3.1, the space-time gradients d/dx^ 
become the four-momenta q^, beginning with Eq. (|65|). In the long-wavelength limit, it is generally 
required that q, which probes the spatial inhomogeneity of the (stationary) system, be small compared 
to the relevant momentum (gradient) scales, i.e. the temperatures T e ,T u . Otherwise, the response 
functions, see Eq. (|80|), would inherit neglected higher order terms in q^d£ , which correspond to going 
from Eq. © to Eq. (g§). 

A truly microscopic transport calculation of the neutrino distribution, as they are released from 
the neutrino sphere, is an interesting topic for future work |l], g]. We wish to conclude by illustrating 
the modifications resulting from a more realistic smooth cut-off neutrino momentum distribution. 

For example, we consider the azimuthal angle integral which contributes the singular term oc 
(q + q ) I ( z q — q°) to the neutrino response function M^ 00 , Eq. @, which in turn is essential for the 
longitudinal Type I pharon originating in Eq. ( |110[ ). Following the radial momentum and polar angle 
integrations, one encounters the integral: 

/ = q 2 J 1 dz (q° + ie-\q)z)- 2 F(z) , (134) 

where we replaced the previous sharp cut-off, i.e. @(z — z m ), by the smooth function F, 

F(z) = iV^T 1 (arctan[\/a(z — z m )\ + arctan[-v/a(l + z m )\) , (135) 
where z m = cos9 max , and we have introduced the convenient abbreviations: 

N± = arctan[y / a(l — z m )] ± arctan[\/a(l + z m )] . (136) 



Thus, we have F(—l) = 0, F(l) = 1 and the Lorentzian derivative, 



Collective Modes in Neutrino 'Beam' Electron-Positron Plasma Interactions 



23 



The parameter a > determines the steepness of the sigmoid cut-off, which ultimately should be 
related to the emission characteristics of the neutrino sphere [^. 

Then, after a partial integration, we obtain an integral with three complex simple poles, 

I -i±M + A[ i dz m {m 

\q\ + WJ-i z-{ q » + l e)m ' 
which can be solved analytically. The final result is: 

Thus the previous result is recovered in the sharp cut-off limit. 

We observe that the previous resonance pole at q° = Qu = \q\z m is moved symmetrically off 
the real axis. However, closer inspection shows that the above result is not singular any more as 
Q° ~> \q\( z m i i/y/a). We recalculated the complete M^) 00 -component of the neutrino response 
function with the smooth cut-off function F(z) and found a corresponding result. 

Therefore, in accordance with intuitive expectation, a physically more appropriate moderate 
smoothing of the angular dependence of the neutrino momentum distribution turns the pharon res- 
onance poles into resonances with a finite width, which decreases with increasing sharpness of the 
cut-off (i.e. as a — ► oo). Most likely, this implies that realistic growth rates for the unstable modes 
will be lower than estimated optimistically in Section 3.2.3 . 



4 Conclusions 

In this work, we studied the collective modes which were earlier conjectured to produce anomalously 
large instabilities in the course of the interaction of the neutrino 'beam' with the plasma sphere in 
Type II supernovae ||. 

For this purpose, we derived the semiclassical transport theory based on the Dirac field equations 
for neutrinos and electrons, which are coupled according to the Standard Model. Our results also allow 
for the handling of situations with strong spin-polarizing magnetic fields, which we did not consider 
here. We derived a related linear response theory and applied it in a detailed supernova scenario, 
adapted from Bethe's review ||. 

We studied, in particular, the modifications of the usual transverse and longitudinal plasmons of 
an electron-positron plasma, which are caused by a high-power beam-like flux of neutrinos, such as the 
almost radially outward streaming neutrinos, which are released from the neutrino sphere surrounding 
the supernova core. In the collisionless approximation, we found only a very weak perturbative effect 
and no induced imaginary part of the dispersion relation. This is due to the suppression of all neutrino 
effects by the small dimensionless effective coupling constant, 

2 _ 1 n2 2 2^ r n-22 

a = ^G F m w m D 10 

under typical Type II supernova conditions. 

However, we found new types of growing, as well as damped collective modes, the pharons, which 
are characterized by an essentially linear dispersion relation tu{q)/q ~ const in the long- wavelength 
limit. Their characteristic properties depend on the relative orientation of the neutrino beam, collective 
mode propagation and electric current fluctuation directions (the Cases I and II are studied in 
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Sections 3.2.2 and 3.2.3). They partially overcome the discussed suppression, since a resonance pole 
arises in the dispersion equations in the region with u>(q) < q. In this region, one finds in an ordinary 
electron-positron plasma strongly Landau damped modes or, in the ultrarelativistic limit, no (plasmon) 
modes at all. 

As we estimated roughly in Section 3.2.3, although they partially overcome the discussed suppres- 
sion, the Type I pharon growth rates are still about four, and likely more, orders of magnitude too small 
to make an impact on supernova evolution. The Type II pharon growth rate is even more suppressed; 
it is proportional to Gp, like a purely weak interaction cross section. Clearly, we have seen that for 
the electro-weak interaction to be relevant, the effective particle densities have to be sufficiently high, 
since one has to overcome the suppression expressed by the smallness of a 2 , given above, or its equiv- 
alent for more massive particles. We remark that changing the neutrino and electron temperatures, 
as compared to the typical supernova case discussed after Eq. (|122j ), for example, while keeping the 
'wavelength' rri£)/q fixed, the Type I pharon damping constants grow and the corresponding e- folding 
lengths drop oc T~ 2 T~ 3 



e 



At this point it is worthwhile remarking that more realistic calculations also have to take the 
collisonal damping into account. This can be incorporated in our approach in the relaxation time 
approximation in future applications ||] . 

Finally, we point out that pharon type modes should occur in other situations with a current- 
current type interaction under two-stream conditions. The anisotropic momentum distribution char- 
acteristic of a 'beam' with limited (momentum) opening angle appears to cause the resonance effect 
exciting these modes by impact on an isotropic plasma. - We did not study here a spatially limited 
beam or jet, which causes quite different 'hydrodynamic' instabilities. 

Pharon modes may perhaps be fed effectively by the still unknown central engine powering gamma 
ray bursts, for example, see Refs. [21]. There is growing evidence for even truly jet-like processes in 
the gamma ray burst phenomenon. Furthermore, allowing for other than electro-weak interactions, 
such modes possibly come into play in the ultimate evaporation of primordial black holes 22]. 

To summarize, we conclude that the intrinsic weakness of the neutrino caused collective effects, 
related to the large asymmetry of the electromagnetic and weak coupling strengths, makes it rather 
unlikely that they play a role in the neutrino energy deposition in the supernova plasma sphere. 
However, the new pharon type instabilities may be quite relevant in two-stream situations occuring 
in other astrophysical systems. 
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Appendix 

Here we consider in more detail the structure of the neutrino-antineutrino spinor Wigner function. 
The general results of Eqs. ©-(0) can be further specialized for the case of the Standard Model 
^L^R-system in the massless limit. - We follow the notation of Ref. |^3] in this Appendix. 

In the main part of the paper we use the Dirac-Pauli representation of the 7-matrices, which are 
defined by the anticommutation relations {7^,7^} = 2<7^ , i.e. explicitly: 




(140) 
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where all entries are 2 x 2-matrices themselves; in particular, a l ,i = 1,2,3 denote the standard Pauli 
matrices. We listed also the chirality operator 7 5 = 2 7 °7 1 7 2 7 3 , which anticommutes with all 7^. In 
the following we will conveniently begin with the chiral Weyl representation: 

•*=(; i) • 0) ■• M - * 1 ;) • <-> 

The subscripts D, W presently serve to distinguish the representations, which are related by: 

7& = Uwl^ w , U w = ^{l-^ wl ° w ) (142) 

i.e. a simple unitary transformation. 

In the Weyl representation the stationary Dirac equation separates into a pair of two-component 
Weyl equations: 

Ex(p> = -o ■ p x(p) , (143) 
Ecj)(p) = +d-p^p{p) , (144) 

the first of which describes the physical v\ ^-system, as we shall see, while the second one is discarded 
in the Standard Model. Using {cr*, cr?} = 25^ , one verifies that E 2 = \p\ 2 , in both cases. Concentrating 
on Eq. ( |143| ) , we consider first the positive energy solution with four- momentum = [E > 0,p), which 
obeys (p = p/\p\): 

°-PX+ = ~X+ • (145) 

It thus describes the negative helicity (left-handed) neutrino, vl. Conversely, for the negative energy 
solution with p^_ = [E < 0, — p), we obtain: 

ff ■ (-p) x- = +X- , (146) 

thus describing a negative four-momentum positive helicity (right-handed) neutrino, which is the 
positive four- momentum right-handed antineutrino, vr. 

We conclude that in the Weyl representation the physical z^i/R-spinor ipM is represented by: 

^=[ X = a % + , a - X -) , (147) 

where x ls written as a suitably normalized superposition of the two-spinors x±i e -g-) considering plane 
wave states [23]. Next, we calculate the corresponding four-spinor in the Dirac-Pauli representation: 

# = ^ } = ^(_ X X ) , (148) 

where we applied the unitary transformation U defined in Eqs. ( |142| ) . 

We proceed to calculate the bilinear covariants which are needed in Eqs. (|7|)-(|T2|). Using the explicit 
form of the 7-matrices in the Dirac-Pauli representation, Eqs. (|140| ), we obtain: 

^("tyM = = ^M 7 5 ^ (l/) , (149) 

with \jj( u ) = ip( u ^j° = (x^ X^)/\/2 ■ dropping the subscript D from now on. Furthermore: 

^Hy^H = - x \j x = ^H^y^H . (151) 
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Finally, using o~^ u = i^y^ ,j u ]/2 , we also obtain: 

^Mff^^M = . (152) 

Summarizing, only the vector and axial vector (two-point) densities are nonzero for the vl^r- 
system and, in fact, they are equal. This yields for the neutrino Wigner function components: 
V/i (x,p) = (x,p), while all other components vanish, see Eqs. (f7|)-(|T2"l). 

Presently we made use only of the algebraic properties of the neutrino spinors. It did not enter that 



the amplitudes a± from Eq. ( 147 ) actually are to be considered as creation/annihilation operators. 
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